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We present our QCD analysis of the proton structure function F^(a;,Q 2 ) to determine the par- 
ton distributions at the next-to-leading order (NLO). The heavy quark contributions to F^{x,Q 2 ), 
with i = c, b have been included in the framework of the 'fixed flavour number scheme' (FFNS). 
The results obtained in the FFNS are compared with available results such as the general-mass 
variable-flavour-number scheme (GM-VFNS) and other techniques used in global QCD fits of par- 
ton distribution functions. In the present QCD analysis, we use a wide range of the inclusive 
neutral-current deep-inelastic-scattering (NC DIS) data, including the most recent data for charm 
F£, bottom F%, longitudinal Fl structure functions and also the reduced DIS cross sections o^r NC 
from HERA experiments. The most recent HERMES data for proton and deuteron structure func- 
tions are also added. We take into account ZEUS neutral current e ± p DIS inclusive jet cross section 
data from HERA together with the recent Tevatron Run-II inclusive jet cross section data from 
CDF and D0. The impact of these recent DIS data on the PDFs extracted from the global fits are 
studied. We present two families of PDFs, KKT12 and KKT12C, without and with HERA 'combined' 
data sets on e ± p DIS. We find these are in good agreement with the available theoretical models. 

PACS numbers: 13.60.Hb, 12.39.-x, 14.65.Bt 



I. INTRODUCTION 

Lepton nucleon deep inelastic scattering (DIS) pro- 
cesses play a very important role for our understanding of 
quantum chromodynamics (QCD) and nucleon structure 
[l| . The DIS process of leptons off nucleons is an instru- 
mental tool for high precision measurements of the quark 
and gluon content of the nucleons. DIS experiments pro- 
vide valuable information from lepton-nucleon scattering 
cross sections to find out about the structure functions of 
targets. Note that these structure functions are related 
to parton distribution functions (PDFs) directly. This 
process at large momentum transfer provides one of the 
cleanest possibilities to test the predictions of QCD and 
allows us to measure the high-precision PDFs of the nu- 
cleons together with the strong coupling constant. DIS 
process also demonstrates the internal structure of the 
nucleon in high-energy scattering at the CERN Large 
Hadron Collider, Fcrmilab Tevatron collider, and in other 
hard scattering experiments. The advent and develop- 
ment of the electron-proton colliders, especially HERA, 
has resulted in significant progress in the recent years on 
the precise understanding of the structure of the nucleon, 
especially in different kinematic regions of momentum 
fractions x of proton carried by the struck parton. The 
most precise inclusive structure function measurements 
which has been provided by HERA collider, has brought 
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remarkable improvements to our understanding on pro- 
ton structure and can be used to extract high precision 
PDFs, which is important for any process which involves 
colliding hadrons. 

Recently it is quite common for several theoretical 
groups to make a global fit on available experimental 
data to extract parton distribution functions which are 
utterly useful in colliding hadrons, specially LHC, stud- 
ies. Nowadays thousands of published data points from 
various experiments arc available to extract more precise 
PDFs and strong coupling constants in higher perturba- 
tive orders. Also for more accuracy, the QCD fits can be 
performed with leading-order (LO), next-to- leading order 
(NLO) and even next-to-next-to-leading order (NNLO) 
QCD approximation. In addition, it is also significant to 
find out about PDFs uncertainties for experimentalist. 

In recent years, a number of theoretical groups 
providing global QCD fits (CTEQ/CT10 @, Q and 
MRST/MSTW08 0-Q) have developed a prescription 
for determining uncertainties, and have now been joined 
by several other theoretical gro ups (NNPDF 043; 
ABKM10 [T(J, GJR08/JR09 [Hill). HERA experimen- 
tal collaborations have combined their DIS data and pub- 
lished parton distribution functions. (HERAPDF), with 
uncertainty error bands as well [HI, [l4| • 

Nowadays, new PDFs incessantly become more accu- 
rate by considering up-to-date theoretical developments 
and data from hadron colliders experiments, also such 
evolvement provide more reliable calculation of uncer- 
tainties corresponding to different theoretical and exper- 
imental inputs. The correct treatment of heavy quark 
flavours in global QCD analysis of PDFs is essential 
for precision measurements at hadron colliders such as 
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HERA, since it contributes up to 30% to the DIS inclu- 
sive cross section at small x. This is very important for 
the reliability of parton distributions in the kinematic re- 
gion corresponding to the foreseen experiments at LHC 
since DIS remains an incomparable source of information 
about parton distributions at small x, despite the effect 
of collider data on the PDFs. Recent reported results 
0, H5M2H show that the cross section of Higgs, W and Z 
production at the LHC are sensitive to detailed aspects 
of PDFs that depend on the heavy quark mass effects; 
and standard model processes as well as beyond stan- 
dard models considerably depend on improved knowledge 
of the c-quark parton density, in addition to the light 
flavours density [22[ • Heavy quark flavours production in 
DIS is calculable in QCD and provides significant infor- 
mation on the gluonic content of the proton. Nowadays, 
heavy quarks arc produced in several experiments of high 
energy physics and the production and decay properties 
of the heavy quarks are extremely interesting. A num- 
ber of theoretical groups, such as CTEQ collaboration, 
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CT10 [3j, MSTW08 [6|, ABKM10 [Io| and GJR08 
have produced publicly available PDFs using different 
data sets. We also studied the diffractive process to pa- 
rameterize the diffractive parton distribution functions 
(DPDFs) taking into account heavy quarks contributions 
using recent experimental data [23j. 



In the present paper we will extract the PDFs using 
QCD analysis of a large number of available and up-to- 
date experimental data in the frame work of FFNS by 
introducing a new kind of parametrization. This paper 
deals with the construction of a FFNS and the compari- 
son of their predictions with the most recent experimen- 
tal data for various region of x and Q 2 . Since we have 
recent experimental data for the reduced cross section, 
heavy quarks and longitudinal structure functions from 
the HERA collaboration, and new fixed target data for 
proton and deuteron structure functions from HERMES, 
we have enough motivation to extract all of parton distri- 
bution functions using new sets of parametrization with 
take into account the heavy quark contributions. 



The organization of the paper is as follows. In Sec. [H] 
we outline the theoretical formalism which describes the 
DIS structure functions. Our parametrization and heavy 
quark contributions to the proton structure function and 
the formulation of FFN scheme are also performed in this 
Section. In Sec. IIIIl we discuss the global parton analysis 
and then present the results of our NLO PDF fit to the 
DIS world data to determine the PDF parameters and 
a s (M§). Discussion of fit results are given in Sec. IIVI 
This section also contains a detailed comparison of the 
KKT12 and KKT12C QCD fits with other recent PDF sets 
and with the recent HERA data in various x and Q 2 
regions. The results related to the longitudinal structure 
function Fl(x,Q 2 ) is presented in Sec. |Vj Section IVIII 
contains our summary and conclusions. 



II. THEORETICAL BACKGROUND OF THE 
QCD ANALYSIS 

In this section, we first give a brief overview of the stan- 
dard theoretical formalism that we used in the present 
QCD analysis. A more detailed description is given later 
in separate sections. 



A. Theoretical input 

The QCD predictions for the structure functions are 
obtained by solving the well known DGLAP evolution 
equations at NLO in the modified minimal subtrac- 
tion (MS) scheme. We choose the factorization and 
renormalization scales to be Q 2 . We work in the MS 
scheme, where the deep-inelastic nucleon structure func- 
tion F2(x, Q 2 ), can be written as a convolution of the par- 
ton distribution functions together the coefficient func- 
tions jf|. The scale dependence of the PDFs is given by 
the well-known DGLAP evolution equations which yield 
the PDFs at all values of Q 2 if they are provided as func- 
tions of x at the input scale of Q 2 ,. 

A number of methods to solve the evolution equations 
have been proposed, including direct x-space method, 
Mellin-transform methods and orthogonal polynomials 
methods. More detailed about the later method are given 
inRefs. . 

In the present paper wc have performed all Q 2 - 
cvolutions in Mcllin A-momcnt space for the NLO evo- 
lution. In this regard, the program QCD-PEGASUS [36| 
is used in order to perform all Q 2 -evolutions of parton 
distributions. The input for the evolution equations and 
the strong coupling constant a s {Q 2 ) ), at a reference input 
scale taken to be Q 2 , = 2 GeV 2 , must be extracted from 
a global QCD analysis of DIS world data. 

Our present standard NLO QCD analysis has been per- 
formed within the MS factorization and renormalization 
scheme. We do not consider the heavy quarks (c, 6, t) 
as massless partons within the nucleon. In this case, 
the number of active (light) flavours rif appearing in the 
Wilson coefficients and the corresponding splitting func- 
tions to be fix at rif = 3. Furthermore, the evaluation 
of the running strong coupling a s {Q 2 ) needs to match 

a ( s ni) aXQ = m h (m c = 1.41 GeV, m b = 4.50 GeV, m t = 

175 GeV), i.e. a<? f \m h ) = ai nf ~ 1] \m h ). 



B. Parametrization 

For the QCD fit in the present analysis we use the fol- 
lowing independent functional form of the parametriza- 
tion of the parton distribution functions at the input scale 
Qq=2 GeV 2 for the valence quark distributions xu v and 
xd v . the anti-quark distributions xS — 2x(u + d + s), 
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xA = x(d — u), and for gluon distribution xg: 

xu v {x,Ql) = A u x a "(l - x)^(l +~/ u x K +rj u x), 

xdvix, Ql) = A d x ad {l-xf"{l + -/ d x Sd +r] d x), 

xA(x,Q 2 ) = A A x aA (l~xf s+ ^{l + lA x s * +ri A x), 

xS(x,Ql) = A s x as (l-x) f3s (l+ ls x Ss + Vs x), 

xg(x, Ql) = A g x a ° (1 - xf' (1 + lg ^ + V g x). (1) 

The flavour structure of the light quark sea at Qq is taken 
to be 

2u(x,Ql) = OAS(x,Q 2 )-A(x,Q 2 ), (2) 
2d(x,Ql) = OAS(x,Q 2 ) + A(x,Q 2 ), (3) 
2s{xMl) = 0.2S(x,Q 2 ), (4) 

and we take a symmetric strange sea s(x, Qq) = s(x, Qq) 
due to the insensitivity of data sets we are using to the 
specific choice of the strange quark distributions. Due 
to the very large amount of data in our global QCD fits 
and the presence of some constrains, we can control the 
parameters in order to allow sufficient flexibility in the 
form of the parton distributions. 

In the present work, the input PDFs listed in Eq.[T]arc 
subject to three constraints from valence quarks number 
sum rules together with the total momentum sum rule 
Q. Not all the parameters in our inputs for the parton 
distributions are free. The parameters A u , A d and A g 
were calculated from the other parameters using men- 
tioned constrains, therefore there are potentially 28 free 
parameters in the PDF fit, including a s (Qo)- 

The mentioned parametrization for the input PDFs in 
our analysis, Eq. [IJ and all free PDF parameters listed 
there allow the fits a large degree of flexibility. In the 
final minimization, we fixed some parameters and only 
13 parameters remained free for all parton flavour in the 
fits. These parameters are allowed to go free in order to 
calculate the covariance matrix. 

According to our parameterization in Eq. [T] and for 
valence quark distributions, xu v and xd v , we use the 
standard functional form. The normalization parame- 
ter A v is obtained from the number of valence quarks 
sum rules. The variation of the valence distributions at 
intermediate values of x will be controlled by the third 
free parameter j v . In our calculation for valence sector, 
we found that the parameter rj v is correlated to a v and 
(3 V more strongly. At very high and even at very small 
x, the uncertainties are constrained by the mentioned 
sum rules. In the region of x =[0.01, 0.75], where the 
valence PDFs are constrained by DIS data, about 75% 
of the total number of valence quarks can be found. We 
should notice that in the reported results of MSTW08 @ 
and GJR08 fill ] for input parametrization, the parameter 
value of 5 V is fixed to 0.5 which we want to deal with this 
parameter as a free parameter. After first minimization, 
in both cases of valence quark distributions, we will get 
the parameter values of 7„, S v and rj v fixed so only two 
parameters remain free, a v and j3 v . 



In order to consider the flavour separation of the sea- 
quark distributions xA = x(d — u) in Eq. [TJ we choose 
As + Pa as a power of (1 — x) for xA, because at high 
x this PDF is becoming very small and we attempt to 
constrain u and d to be positive. We need also to have 
x 5& term instead of the more usual x 05 and x 2 term to 
give sufficient flexibility at medium and large x. In fact 
we consider A a, a a and /3 a as free parameters and 7a, 
<5a and r)A will be fixed at their best-fit values. 

In the sea quark distribution, xS — 2x(u + d + s) in 
Eq. [U there is no sum rule constraint, so in principle 
there are 3 free parameters including As, as and /3s and 
the other 3 parameters will expected to be fixed since the 
data do not constrain these parameters well enough. 

The gluon distribution is a far more difficult case for 
PDFs parameterizations to obtain precise information. 
Since the current data provide little constraint on the 
gluon density, we expect the functional form of gluon 
distributions in Eq. [1] to allow us for better gluon PDF 
determination. The input gluon in the MSTW08 His 
almost the same form of the original MRST fits (37U38T] , 
but MSTW08 added a second term to avoid a negative 
gluon distribution at very small x. We found that the 
parameter S g in our gluon parametrization could help us 
to facilitate this. In this case we need to have a x bs=2 
term instead of the more usual \fx term to give sufficient 
flexibility at low x, and we find stability in the QCD fits. 
On the other hand, one feature of our parametrization 
is that we allow the gluon to have a very general form 
over a wide range of x, specially at medium x, because it 
can induce potentially large effects on the strong coupling 
constant obtained from the global fits and indirectly leads 
to a change in the light partons. In particular, fixing the 
parameters 7 9 , S g and rj g at the best fit values we find that 
a g and /3 g are sufficient for the gluon - one for high x and 
one for low x. Consequently in the final minimization, 
only two parameters a g and (3 g remain free. We note 
that we can determine A g by the total momentum sum 
rule @]. 

As a last point, we note that the single-inclusive jet 
data are very important for the constrains of gluon dis- 
tribution. The constraints on the gluon distributions in 
the domain 0.01 < x < 0.5 are provided by the inclusive 
jet production data from Run II at the Tevatron from 
CDF and D0 while the data from HERA on inclusive jet 
production in DIS constrain the gluon for 0.01 < x < 0.1. 



C. Treatment of heavy flavours 

In this section, we will discuss the fixed flavour number 
(FFN) scheme parton model predictions at high energy 
colliders at LO and NLO of QCD. It is expected that 
perturbative QCD at NLO should give a good descrip- 
tion of heavy quark flavour production in hard scattering 
processes, so the correct treatment of heavy flavours in 
common QCD analysis of PDFs at this order is essential 
for precision measurements at hadron colliders. The LO 
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and NLO heavy quark contributions F°' (x, Q 2 ) to the 
proton structure function are calculated in the FFN ap- 
proach and contribute to the total structure functions as 
Fi(x, Q 2 ) = if ght + if cavy where the if ght refers to the 
common it, d, s (anti) quarks and gluon initiated contri- 
butions, and if cav y = F? + F? with i = 2,L. Charm 
and bottom quark contributions to the total DIS cross 
section at high Q 2 are up to 30% and 3%, respectively, 
and top quark contributions are negligible at present en- 
ergies. As already mentioned at the end of Sec. Ill A[ 
we employ for our analysis the FFN scheme and fix the 
number of active light flavours rif = 3 in all splitting 
(k) 

functions and in the corresponding Wilson coeffi- 

cients. In this factorization scheme only the light quarks 
(u, d, s) are genuine, i.e., massless partons within the 
nucleon, whereas the heavy ones (c, b, t) are not. 

Fixed flavour number approach is the region where the 
hard scale of the DIS process is around the quark mass, 
i- G - Q 2 ^ m h where h = c,b, t. In this case, massive 
quarks are not considered as partons within the proton. 
They are described as final-state particles and so they are 
entirely created in the final-state. Our practical reason 
for extracting partons from the fits to structure function 
data in the fixed three- light-flavour (n / = 3) is to have a 
set of up-to-date partons in order to make consistent com- 
parisons between other parton sets available using the 
same theoretical framework of other scheme such as gen- 
eral mass variable-flavour number scheme (GM-VFNS) 
[H, Q and FFNS 11 1. In addition for many exclusive 
or semi-inclusive processes the theoretical predictions for 
the hard scattering cross sections are only available in 
the FFN scheme. The different heavy quarks treatment 
which nearly every group does, may led to different re- 
sults for the extracted parton distribution functions. 



with the non-singlet contribution for three active (light) 
flavours are given by 



1 + 



r (o) r (i) 
L , 1 



(x,Q 2 ), (6) 



where C^ q {z) = 5(1 — z) and C^NS * s *hc common NLO 



coefficient function [4J-|44| ■ The NLO (^-evolution of 
the flavour non-singlet combinations q% = u + u + d+d— 
2(s + s) = u v +d v + 2u+2d— 4s and gy = u + u— (d+d) = 
u v — d v is related to the LO (1-loop) and NLO (2-loop) 
splitting functions P^s and i^g respectively [42|, |43[ . 
The flavour singlet contribution in Eq. [5] reads 



-F 2t s(x,i 
x 



■>-§{ 



a 



(0) 
2,9 



« C. ■ 



with S(cc,Q 2 



/}(^,Q 2 ), 



(7) 



a 



(i) 

2,9 



a 



(i) 

2,NS 
.(1 



J + d v + 2u + 2d + 2s, 
and the additional NLO gluonic coefficient 

function C^g can be found in Refs. [HI, [44[ for example. 
Notice that we consider sea breaking effects (u ^ d) since 
the HERA data is used, and thus our analysis are sensi- 
tive to these corrections. Since these data sets are also 
insensitive to the specific choice of the strange quark dis- 
tributions, we consider a symmetric strange sea, s = s. 
We also use s = s=^(u + c[) where in practice k is a 
constant fixed to k = 0.4 — 0.5. The heavy flavour contri- 
butions F£ ' b (x, Q 2 ) that have been added in the present 
analysis will be discussed in the next section. 



E. Heavy flavour contributions 



D. Fixed flavour number scheme (FFNS) and F2 
structure function 



The heavy flavour contributions F2' b (x,Q 2 ) which 
have been added in the present analysis are taken as in 
Refs. (4lL |45|-|50| . The structure functions in fixed flavour 
factorization scheme are written by 



The NLO QCD analysis presented in this paper is per- 
formed in the modified minimal subtraction (MS) factor- 
ization and renormalization scheme. Heavy quarks (c, 
b, t) are not considered as massless partons within the 
nucleon and all their effects are contained in the pcr- 
turbativc coefficient functions. This defines the so-called 
"fixed flavour number scheme" (FFNS), which is fully 
predictive in the heavy quark sector and incorporates the 
correct threshold behavior. In the common MS factoriza- 
tion scheme, the relevant total proton structure function 
F%(x, Q 2 ) as extracted fr om the DIS _ep process can be, 
up to NLO, written as 



irom tne Dib ep 



F 2 (x,Q 2 ) = F 2 \ s (x,Q 2 )+F 2 , s (x,Q 2 ) 



2,NS 

Ft b \x,Q 2 ,ml b ). 



(5) 



f£(x,q 2 ) 



E 



Cr™ S ' nf (Q 2 /ml)®f? f (Q 2 ), (8) 



where f™ f (Q 2 ) represents only to the light partons (i = 
u, d, ...), rrih refers to the mass of the heavy quarks, with 
h = c , b and nj denotes the number of active flavour. 

The FFNS approach which we are using in the present 
analysis, contains all the heavy quark mass dependent 
contributions. In leading order of QCD, 0(a s ), the well 
known FFN scheme heavy quark flavour contributions 
C^™ S ' nf {Q 2 /m 2 ,) in ep collisions are dominated by the 



photon gluon fusion process (BGF), 7*3 — > hh [5JH55J. 
The FFNS heavy flavour coefficient functions for neutral- 
current structure functions containing ln(Q /m?) terms 
arc well known up to NLO. Due to the complexity of the 
heavy flavour coefficients function caused by the non-zero 
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mass of heavy quarks, nth, there arc no published com- 
plete analytic expressions for all the coefficient functions 
and the results are available in table form at this order 
[49l | . In particular we use the code provided in Ref. [4!| , 
which combines known analytic expressions together with 
grids for the more complicated coefficient functions; they 
are represented in the MS factorization scheme. Up to 
0(a 2 ), the relevant heavy quark contributions to DIS 

are given by (4!| [Ho| 



structure functions F 2 ^ vy 



F^ =2L (x,Q 2 ,m 2 ) 



Q 2 a 2 



Q 2 a s 



47r 2 m 2 
dz 



dz 



x ,.2\Jfi) 



z 



2 



i=q,q 



(9) 



where the upper boundary on the integration is given 
by Zmax = Q 2 /(Q 2 + 4m 2 ) and sum runs over all 
the light quark and antiquark flavours. Furthermore, 
fi(x,/i 2 )(i = q,q) denote the light parton densities in 
the proton and fi stands for the mass factorization scale, 
which normally has been put equal to the renormal- 
ization scale, is the charge of the heavy quarks, 
with h = c,b. The coefficient functions, represented by 

c%,0,Cj}(r),0,(i = 9,Q,q;l = 0,l;j = 2,L) and 
by df}( V ,0,$Uv,0, (i = = 0,1; j - 2,L) are 

calculated in [50[ and they are represented in the MS 
scheme. The LO 0(a s ) contributions to F 2L , due to 
the subprocess j*g hh with h = c, 6, have been sum- 
marized in [55[, and the NLO 0(a 2 ) ones are given in 

mm. 

Different groups use various values for the charm and 
bottom quark masses varying from 1.3 GeV to 1.65 GeV 
and from 4.3 GeV to 5 GeV, respectively. This variation 
will change the extracted parton distribution functions. 
A detailed discussion on mass dependence of parton dis- 
tribution functions can be found in [jj d, [H, [57| ■ 

In the present analysis we fix the heavy quark masses 
at m c = 1.41 GeV and mi, = 4.50 GeV, this differs from 
the MRST and MSTW08 default values of m c = 1.43 
GeV, m b = 4.30 GeV and m c = 1.40 GeV, m b = 4.75 
GeV, respectively. The results of Ref. 0] indicated that 
if we allow the mass of charm quark as a free parameter 
in the global fits, one can find the best-fit value to be 
m c = 1.39 GeV at NLO. 



III. GLOBAL PARTON ANALYSES 
A. Overview of data sets 

PDFs are usually determined in a global QCD analy- 
ses of a wide variety of DIS data, both from HERA and 



fixed-target experiments as well as v(v)N xF% data from 
CHORUS, NuTeV and data for the longitudinal structure 
function, Fl(x, Q 2 ). In the recent years, several new pre- 
cise data sets became available which we include in our 
present global QCD fit. In the present analysis we fit to 
the measured neutral current deep inelastic e ± p scatter- 
ing cross section values that is defined by three structure 
functions, F 2 , Fl and xF 3 as fUfllj 



VNC 



(x,Q 2 ) 



F? c (x,Q 2 )- 



F^(x,Q 2 )T^F, 
1 + 



xc 



(10) 
Q 2 /(xs) 
1±(1- 



where < y < 1 is the process inelasticity y 
and s is the ep centre-of-mass energy and Y± 

v 2 )- 

Statistically, most significant data that we use in our 
QCD analysis, which we call "KKT12" fit, are the HERA 
measurements of the DIS 'reduced' cross-section. We use 
the small and large- a; from HI and ZEUS <J^(x, Q 2 ) data 
for Q 2 > 2 GeV 2 [HH . 

In our QCD fits, we use the most recent HI and ZEUS 
combined measurement for inclusive e^p scattering cross 
sections [l4| instead of all the above mentioned data. To 
distinguish, we call this "KKT12C" when we used the re- 
cent HI and ZEUS combined measurement. 

This combined data set contains complete informa- 
tion on DIS cross sections published by the HI and 
ZEUS collaborations. The kinematic range of the neu- 
tral current data is 0.045 GeV 2 < Q 2 <30000GeV 2 and 
6xl0 -7 < x <0.65, for values of inelasticity y between 
0.005 and 0.95. The statistical uncertainty reduces in 
this measurement. 

The most recent HERA structure function data which 
we use in our analysis are the HI measurements of the 
e^p DIS 'reduced' cross-section (69|. The kinematic 
range of the following measurement covers 1.5 GeV 2 < 
Q 2 <120GcV 2 and 2.9xl0" 5 < x < 0.01. By using 
these data, we can improve our PDFs, especially at small 
x value. 

In addition, we have used the fixed tar get F% data of 
NMC [71, BCDMS [zj, E665 [z2, SLAC [23 and HER- 
MES [zj all subject to the standard cuts Q 2 > 2 GeV 2 
and W 2 = Q 2 {\ - 1) + mp>. 12.5 GeV 2 . Furthermore the 
fixed target data for F^j|3, [ll III III and the structure 
function ratio F 2 /F 2 [76l | are both subject to the same 
standard cuts have been added to the present analysis. 
These data on ff , F^ and F 2 d /F| help us to control the 
behavior of PDFs in large x, which dominated by valence 
quarks [l[. 

The HI and ZEUS measurements on the heavy flavour 
contribution to structure functions F 2 , F 2 and F 2 /F 2 
of [77l - l88| have also been included in the both sets of 
NLO fits. The Drell-Yan dimuon pair production data of 
the E866/NuSea (fixed target) experiment [H, |90[ have 
been used. 

The direct HERA measurements of the longitudinal 
structure function, F^ [go|, HH, HH, |9l| which is a direct 
constraint on the small x region of gluon distribution, 
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- 30000.0 


147 [61] 


1.5% 


0.9999 


HI Q 2 99-00 e+p NC 




o.Orj-D 


n no 
- O.Uz 


2.0 - 


12.0 


OO iDZl 


U.o/o 


0.9994 


HI Q 2 00 e+p NC 






n 1 
- U.l 


12.0 - 


- 150.0 


GO TfiQl 

yy [0^1 


U.o/o 


U.yyo ( 


ZEUS SVX 95 e+p NC 




5.4E-5 


- 0.0019 


2.50 - 


- 17.0 


30 [64] 


1 crO/ 
1.070 


1 


ZEUS 96-97 e+p NC 




6.32E-5 


- 0.65 


2.70 - 


- 30000.0 


242 [60J 


00/ 
27o 


0.9995 


ZEUS 06-07 e+p NC 




5.0 E-4 


- 0.007 


20.0 - 


- 130.0 


162 [66J 


1-70/ 

2.7% 


1.0012 


ZEUS 98-99 e"p NC 




0.050 - 


0.65 


200.0 


- 30000.0 


92 [67] 


1.8% 


0.9996 


ZEUS 99-00 e+p NC 




0.050 - 


0.65 


200.0 


- 30000.0 


90 [68] 


2.5% 


0.9968 


HI 03-07 e+p NC 




2.9E-5 


- 0.01 


2.0 - 


120.0 


134 [69] 


4 % 


0.9997 


NMC fip F 2 




0.0045 - 


- 0.5 


2.50 - 


- 65.0 


126 [70] 


2 % 


1.0023 


NMC fid F 2 




0.0045 - 


- 0.5 


2.50 - 


- 65.0 


126 [70] 


2 % 


1.0023 


NMC fJ,n/fxp 




0.008 - 


0.675 


2.23 - 


- 99.03 


156 [76] 


0.15% 


1 


BCDMS up F 2 




0.07 - 0.75 


7.5 - 


230.0 


167 [71] 


3% 


0.9900 


BCDMS lid F 2 




0.07 - 0.75 


8.75 - 


- 230.0 


155 [75] 


3% 


0.9900 


E665 up F 2 




0.0037 - 


- 0.38726 


2.046 


- 64.269 


53 [72] 


1.85% 


1.0012 


E665 nd F 2 




0.0037 - 


- 0.38726 


2.046 


- 64.269 


53 [72] 


1.85% 


1.0012 


SLAC ep F 2 




0.007 - 


0.65 


2.01 - 


- 22.21 


53 [73] 


2% 


1.0128 


HERMES ep F 2 




0.006 - 


0.66 


2.0 - 


12.0 


39 [74] 


7.55% 


1.0167 


HERMES ed F 2 




0.006 - 


0.66 


2.0 - 


12.0 


39 [74] 


7.55% 


1.0167 


H1 en F£ 




1.3E-4 


- 0.02 


6.50 - 


- 60.0 


25 [77] 


7.5% 


1 


HI ep F| 




0.005 - 


0.032 


200.0 


- 650.0 


4 [78] 


1.5% 


1 


HI ep F 2 C 




2.0E-4 - 


- 0.05 


5.0 - 


2000.0 


29 [79] 


1.5% 


1 


HI ep F 2 C 




1.97E-4 


- 0.05 


12.0 - 


- 60.0 


6 [80] 


1.5% 


1 


HI ep F 2 C 




1.3E-4 - 


- 0.00316 


3.50 - 


- 60.0 


10 [81] 


1.5% 


1 


HI ep F 2 C 




8.0E-4 


- 0.008 


12.0 - 


- 45.0 


9 [82] 


1.5% 


1 


ZEUS ep F 2 C 




1.3E-4 


- 0.00676 


4.20 - 


- 111.8 


5 [83] 


1.8% 


1 


ZEUS ep F 2 C 




1.3E-4- 


0.02 


4.0 - 


130.0 


18 [84] 


1.65% 


1 


ZEUS ep F 2 C 




3.0E-5 - 


- 0.03 


2.0 - 


500.0 


31 [85J 


2.2% 


1 


ZEUS ep F? 




8.0E-5 - 


- 0.03 


30.0 - 


- 1000.0 


8 [881 


1 5% 


I 


HI ep F| 




0.005 - 


0.032 


200.0 


- 650.0 


4 [781 


1.5% 


1 


HI ep F 2 b 




2.0E-4 - 


- 0.05 


5.0 - 


2000.0 


12 [79J 


1.5% 


1 


HI ep F 2 b 




1.97E-4 


- 0.05 


12.0 - 


- 60.0 


U OUJ 


1.0/0 


1 

_L 


ZEUS ep F| 




2.0E-4 


- 0.013 


12.0 - 


- 600.0 


9 [86J 


2 % 


1 


ZEUS ep F| 




1.3E-4 


- 0.013 


3.0 - 


450.0 


11 [871 


2 % 


1 


ZEUS ep F 2 6 




8.0E-5 - 


- 0.03 


30.0 - 


- 1000.0 


8 [88] 


1.5% 


1 


HI ep F2/F2 




8.0E-4 


- 0.008 


12.0 - 


- 45.0 


9 [82 


1.5% 


1 


FNAL E866/NuSea 


0.026 - 


0.315 


54 (F 


ixed) 


on Ton nnl 

30 [89, 90] 


0.6% 


1 


Hl/ZEUS F L 


4.27E-5 


- 0.0049 


2.0 - 


110.0 


127 160, 66, 69, 91] 




1 


CHORUS vN xF 3 




0.02 - 0.65 


2.052 


- 81.55 


oO 1921 


010/ 


1 nnoo 
1.UU2J 


NuTeV z/AT 1F3 




0.015 - 


0.75 


3.162 


- 125.89 


64 [93] 


2.1% 


1.0023 


D0 I pp incl. jets 












90 [99] 


6 % 


1.0021 


CDF II pp incl. jets 












76 [100] 


5.8% 


1.0023 


D0 II pp incl. jets 












110 [101J 


6.1% 


1.0221 


ZEUS 96-97 e+p incl. jets 










30 [97] 


2 % 


0.9988 


ZEUS 98-00 e+p incl. jets 










30 [98] 


2.5% 


0.9980 


CDF II pp incl. jets 












20 [1021 


6 % 


0.9972 


All data sets 






3279 







TABLE I: Data sets fitted in our NLO QCD analysis. The fitted normalisations ATn of the data sets included in the global 
QCD fit, together with the total normalisation uncertainty, AJ\f n , for each data set n are also shown in the table. The details 
of corrections to data and the kinematic cuts applied are contained in the text. 
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have also been included in the present analysis. And fi- 
nally, we include the v{v)N xF 3 data from CHORUS H3 
and NuTeV [93| with the same cuts in the present anal- 
ysis. The importance of xF 3 data is that the small dif- 
ference between NC e + p and e~p cross sections data at 
high Q 2 , which determines 2F3 as seen in Eq. [101 contain 
the valence quark distributions. The reason is that the 
structure function xF 3 is related to valence quark contri- 
butions and it can provide direct tests of the total valence 
quark distribution. 

We now indicate that we use the deuteron structure 
function data in our QCD analysis also. Experimen- 
tal data on the deuteron structure function F% and the 
structure function ratio F^/F? have to be corrected for 
nuclear shadowing effects [94I496I ] . Since the neutron 
DIS data come from the experiments with nuclear targets 
such as deuteron, these effects have to be taken into ac- 
count. Taking the shadowing into account, the deuteron 
structure function F^{x^Q 2 ) is related in the following 
way to F%{x, Q 2 ), F^ix^Q 2 ) and to the shadowing con- 
tribution 5Fg(x,Q 2 ), 



B. The method of the minimization and error 
calculation 

1. Minimization of x 2 

The global goodness-of-fit procedure follows the usual 
chi- squared method with X global (p) defined as [1, [H], HH, 
IT07] 



Xelobal (P) 




F*(x,Q 2 



F$(x,Q 2 )+F?(x,Q 2 ) 



,(12) 



where p denotes the set of 13 independent parameters 
in the fit, including a s (Qo)i and 

n data is 

the number of 

data points included, n data = 3279 (KKT12) and 2485 
(KKT12C) for the NLO fit. The errors include system- 
atic and statistical uncertainties, being the total exper- 
imental error evaluated in quadrature. For the i th ex- 
periment, Ff ata , AFf ata , and T" loory denote the data 
value, measurement uncertainty (statistical and system- 



-SF^ix, Q 2 ) . (11) gx\ c combined), and theoretical value for the n th data 



We study the shadowing effect as a nuclear effect in our 
QCD fits and we found that in the range of 0.1< x <0.75 
for the NMC, BCDMS and E665 F^and F 2 d /F| data, 
the shadowing effect are very small [96| . 

The HERA [13, HH and Tevatron [99l-fl03t Run I & II 
data on inclusive jet production in e + p and pp scattering 
have been used in our QCD analysis to extract the new 
sets of parton distribution functions. Take into account 
the jet production data, constrains the the mid- to high- 
x {x « 0.01 — 0.5) gluon density and allows an accurate 
extraction of strong coupling constant a s (A/§) at NLO. 

The fastNLO package [l04j], based on NL0JET++ fiol 
Il06j , presents a method that offers an exact and very fast 
pQCD calculation for a large number of jet data sets. 
This scenario also allows the inclusion of the NLO hard 
cross section corrections to both the Tevatron and HERA 
jet data in the fitting procedure. In the present PDF fits, 
we use the fastNLO code to take full advantage of direct 
sensitivity of the jet production data to the gluon density 
in the proton. 

The total data sets that we use in the present global 
analysis are listed in Table U and are ordered according 
to the type of process. The fitted normalization Af n of 
the data sets included in the global fit, together with the 
total normalization uncertainty, AjV n , for each data set 
n are also shown in this table. There are 3279 and 2485 
data points for KKT12 and KKT12C fits, respectively. We 
note that the data from Refs. [69| to [93| are used in both 
of KKT12 and KKT12C QCD fits. It can be a new feature of 
our analysis that we use a wide range of the most recent 
DIS data sets including different processes, NC, e + p and 
e~p scattering to extract all PDFs. 



point for proton structure function as an example. The 
theory prediction T* eory depends on the input PDF pa- 
rameters p. The KKT12 and KKT12C fits added all uncorrc- 
lated and correlated experimental errors in quadrature in 
the Xgi oba i definition. AJ\f n is the experimental normal- 
ization uncertainty and Af n is an overall normalization 
factor for the data of experiment n. 

In the present analysis, we allow the normalization of 
several data sets to be free at the same time as the PDF 
parameters, and then fixed them in the determination of 
the PDF uncertainties. We also allow a different nor- 
malization for each of the HERA running periods. To 
determine the best fit values at NLO, we minimize the 
Xgiobai w ith respect to 27 free input PDF parameters, 
together with as (Qo) and 23 different data set normal- 
izations, giving a total of 50 free parameters in the fits. 
The minimization of the above \ 2 value to determine the 
best parametrization of th e pa rton distributions is done 
using the program MINUIT jl08l |. Due to the limitation of 
MINUIT package in the practical maximum number of free 
parameters, we fixed some of the data set normalization 
to be 1. For example, in the KKT12 fit, we take all the 
normalization of the HI and ZEUS data on F| c (a;, Q 2 ) to 
be fixed at 1. (See Table [T]) As discussed in Sec. IIII Al 
we use the Hl/ZEUS combined data in our KKTC fit, con- 
sequently the number of free normalization decreases so 
we have this chance to allow the F| c (x, Q 2 ) data to have 
different normalization. 

In the global fits, we obtain 

X 2 3590.589 



NDF 



3266 
2622.118 

2472 



1.098 (KKT12) , 
1.060 (KKT12C) 



(13) 



with the total data sets in Table U and for the parameter 
values listed in Table IIIII Consequently the x 2 fit value 
for KKT12 is almost equal to KKT12C QCD fit at NLO. 



2. Error calculation using Gaussian method 

The one a errors for the parton density xfi(x, Q 2 ) are 
given by the Gaussian err or propa gation method as dis- 
cussed in details in Rcfs. In this method, the 



standard linear errors propagation allows one to calculate 
the error on any quantity F(pi) using the formula 



AF 



dF „ , ,dF 



fri fri dp* d Pi 



"I 1/2 



(14) 



where Ax is the allowed variation in x an d Cjj(p) = 
(H^ 1 )^ is the covariance, or error matrix of the parame- 
ters determined in the QCD analysis at the scale Q\ and 
n = 13 is the number of free input parameters in the 
global fit (see Table IE]). 

Consequently one can calculate the error of any quan- 
tity by using Hessian or covariance matrix. This method 
recently used in our nucleon pola rized and un-polarized 
structure function QCD analysis [ 1 1 lN 1 1 3 1 . The covari- 



ance matrix elements for 13 parameters in the KKT12 fit 
are given in Table HI1 We able to calculate the uncer- 
tainties of PDFs using these covariance matrix elements 
based on the Gaussian method as mentioned in this sec- 
tion. Their value at Q 2 are calculated by the QCD evo- 
lution. 



IV. DISCUSSION OF FIT RESULTS 



Both KKT12 and KKT12C fits show acceptable agree- 
ments with the data sets listed in Table |U In order to 
discuss the fit results, we will concentrate on KKT12 fit. 
Since the outcome of the KKT12 and KKT12C for x 2 is al- 
most similar and for PDFs is slightly different, we will 
show the results of both QCD fits in some figures sepa- 
rately. 



A. Detailed comparison to the HERA and fixed 
target data 

Representative comparisons of our results with the rel- 
evant HERA (HI and ZEUS) data |HM H [g| on 
the structure function of the proton, F%(x, Q 2 ), are pre- 
sented in Figs. Q] and O As shown, the data are well 
described throughout the whole medium- to small-x re- 
gion for Q 2 > 1.5 GeV 2 , and thus perturbative QCD is 
here fully operative. It should be emphasized that the 
perturbatively stable QCD predictions are in agreement 
with all recent high statistics measurements of the Q 2 - 
dcpendancc of F2(x, Q 2 ) in wide range of x. 
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FIG. 1: Comparison of our standard NLO (MS) results for F%(x,Q 2 ) with HERA data EMtfL IftH [65j and Dortmund group, 
GJR08, [llj . To facilitate the graphical presentation we have plotted F^x, Q 2 ) + 0.4 x i(Q ) with i(Q 2 ) indicated in parentheses 
in the figure. 
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FIG. 2: As in previous figure but for large values of Q 2 and larger x. 
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A detailed comparison of our next-to-leading order 
(NLO) (MS) results for F%(x 7 Q 2 ) with a selection of 
HERA data [siM H El and fixed target data of NMC 
0, BCDMS[7l|] and SLAC [z3 are shown in Fig. 
To facilitate the graphical presentation we have plotted 
F^{x, Q 2 )+c with c indicated in parentheses in the figure. 

A comparison of our KKT12C theory predictions for the 
reduced cross section with the Hl/ZEUS combined NC 
e + p and e~p data [l4| are shown in Figs. 2] [5] and [B] 
The error bars in the figures indicate the total experi- 
mental uncertainty. From the figures it is clear that our 
KKT12C fits are in acceptable agreement with the com- 
bined Hl/ZEUS set of reduced DIS cross sections. 

Fig. [7] shows our prediction for reduced cross section, 
Cr,Nc{x, Q 2 ) as a function of x and for different value of 
Q 2 , in comparison with the most recent data from HI 
collaboration 1691. 



B. 



KKT12 and KKT12C PDF sets 



The resulting parameters of the KKT12 and KKT12C fits 
are summarized in Table Hill The x 2 /dof for both PDFs 
fits also shown as well. As mentioned in Sec. IIIB1 the 
values without errors have been fixed after a first mini- 
mization since the data do not constrain these parameters 
well enough. The parameter errors quoted are due to the 
propagation of the statistical and systematic errors in the 
data. 

The KKT12 and KKT12C fitted parton distribution func- 
tions at NLO for all sets of parametrizations and their 
errors at the starting scale Q 2 , = 2 GeV 2 are presented in 
Figs. [5] and |H1 The colored bands show uncertainties on 
these parton distributions. By having the covariance ma- 
trix, it is possible to determine the uncertainties of our 
PDFs for both fits. Although some of the extracted PDF 
parameters of two fits arc a little different but in total we 
found no big differences between the best-fit solutions of 
the two PDF fits. But the error bound for some PDFs of 
two sets are different. 



1. Comparison with other recent PDF sets 

Our best-fit KKT12 PDFs for u-, d-. s-quark and gluon 
distributions xg are shown in Figs. 1101 together with a 
comparison of our results with CTlo||, MSTW08 @, 
ABKM10 [Hj and GJR08 [ll|. 

In Fig. [11] the valence quark distributions xu v {x,Q 2 ) 
and xd v (x, Q 2 ) are shown as a function of x at Q 2 = 100 
GeV 2 in the NLO approximation which have been com- 
ared with the results obtained by CT10 11 , MSTW08 
, ABKM10 [l0|, GJR08 [HI and KT08 H. 
The NLO results for d — u and d/u as a function of 
x in comparison to the results obtained by CT10 3, 
MSTW08 [|, ABKM10 [HJ and GJR08 [TO] have been 
shown in Fig. QH The E866/NuSea results [H, HI. 
scaled to fixed Q 2 =54 GeV 2 , arc shown as the circles 



Parameter 


Fit A - (KKT12) 


Fit B - (KKT12C) 


A u 


0.3346 


0.3753 


Ot u 


0.3271 ± 0.0039 


0.3126 ± 0.0041 


Pu 


3.5858 ± 0.0123 


3.6536 ± 0.0103 


Su 


0.4573 


0.5036 


In 


5.4618 


3.2716 


Vu 


20.1864 


21.2864 


A d 


0.3730 


0.6124 


Ctd 


0.3668 ± 0.0105 


0.4226 ± 0.0092 


Pd 


4.8277 ± 0.1001 


5.1529 ± 0.0654 


S d 


0.6306 


0.9553 


Jd 


4.6633 


-3.7475 


Vd 


7.5520 


14.1380 


A A 


10.0048 ± 0.8417 


15.8096 ± 0.3160 


QA 


1.7030 ± 0.0290 


1.4430 ± 0.0336 


Pa 


9.8938 ± 0.3506 


10.4343 ± 0.4278 


Sa 


0.6625 


0.4972 


7A 


-7.3512 


-7.5151 


VA 


23.7067 


20.4901 


As 


0.3494 ± 0.0019 


0.3506 ± 0.0023 


as 


-0.1740 ± 0.0012 


-0.1775 ± 0.0012 


Ps 


8.2970 ± 0.0953 


8.9898 ± 0.0990 


Ss 


0.5129 


0.8764 


■ys 


0.2434 


0.4549 


T)S 


9.8585 


12.4246 


A g 


6.2512 


7.1383 




nil vo -i- n m m 


U.lo4o ± U.UUol 


P 9 


8.4027 ± 0.1663 


10.0188 ± 0.1494 


la 


0.2600 


1.3264 


Vg 


-1.4998 


-0.9499 


xVdof 


3590.589/3266 = 1.098 


2622.118/2472 = 1.060 



TABLE III: Minimum values of x 2 together with the input 
PDFs parameters at Q% — 2 GeV 2 determined from the two 
different global analysis (Fit A for KKT12 and Fit B for KKT12C 
case) . 



with statistical and systematic uncertainties. The ac- 
curate measurements from the E886 /NuSca experiments 
from 0.02 < x < 0.32 give clear evidence of the d — u 
asymmetry as seen in Fig. 1121 The asymmetry seems 
to reach a maximum at x « 0.2, and becomes small at 
smaller x values. 



2. Impact of the Hl/ZEUS Combined data 

In order to investigate the effect of Hl/ZEUS combined 
e ± p reduced cross section data on the extracted parton 
distributions functions, we show the ratios of the PDFs, 
R^Q 2 ) - AMx,Q 2 )/Mx,Q 2 ) = PDF ' pD -™/< , 
in the central PDF sets of the KKT11C and KKT11 fits 
in Figs. [13] and [H at Q 2 , = 2 GeV 2 and Q 2 = 100 
GeV 2 , respectively. The inclusive jet data from HERA 
and Tevatron Run-II are no t include d. The following re- 
sults can be found in Refs. [HI El- Clearly the effects 
of the Hl/ZEUS combined data sets is on the behavior of 
gluon and valence quark PDFs at x below around 10 _1 
as the Fig. [TBI shows. The KKT11C gluon distribution at 
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FIG. 3: Comparison of our standard NLO (MS) results for F 2 p (x, Q 2 ) with HERA data [5a-(6l|, |&f [65fl and fixed target data of 
NMC [73], BCDMS [zH and SLAC [zl]. To facilitate the graphical presentation we have plotted F%(x, Q 2 ) + c with c indicated 
in parentheses in the figure. 



input scale Q 2 , is dominant for small x, i.e. x < 0.1 while 
the KKT11 u- valance quark PDFs are slightly enhanced at 
very small x, while the (i-valance quark PDFs degraded 
at this region. 

In this region, the KKT11 u, d and s distributions 
are very close to the corresponding KKT11C distributions, 
therefore Afi(x,Q 2 ) — > 0. In addition the sea-quark (u- 



, d- and s-quarks) PDFs are slightly similar at small to 
large x region for these two PDF sets. 

Fig. [T3] shows how these ratios are impacted by the 
DGLAP evolution to higher energy at Q 2 = 100 GeV 2 . 
At small to large x, the ratios for all quarks and the 
features of PDFs are similar to those at Qq = 2 GeV 2 
described above and have been shown in Figures. The 
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FIG. 4: Our results of KKT12C fit for reduced cross section, (j^ Nt ^(x 1 Q ), in comparison with HERA combined NC e^p reduced 
cross section [3] as a function of x for different values of Q 2 bins for 2 GeV 2 < Q 2 < 120 GeV 2 . The error bars indicate the 
total experimental uncertainty. 



same behavior of the u- and rf-valancc quark PDFs at 
Qq = 2 GeV 2 is repeated for higher scale, Q 2 = 100 
GeV 2 . 

The differences observed between the two KKT11C and 
KKT11 PDFs sets using the Hl/ZEUS combined and the 
separate HERA data sets are very small, except for the 



gluon PDFs. At higher scales we have ,niJ — > as 
Fig. fT4l shows. 

The effect of the Hl/ZEUS combined e ± p reduced 
cross section data including the inclusive jet data on ex- 
tracted PDFs from the present QCD analysis have been 
shown in Figs. IT51 and 1161 The results presented in these 
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FIG. 5: Our results of KKT12C fit for reduced cross section, c^nc^' ^ comparison with HERA combined NC e^p reduced 
cross section [l4| as a function of x for different values of Q 2 bins for 150 GeV 2 < Q 2 < 30000 GeV 2 . The error bars indicate 
the total experimental uncertainty. 



figures inclu des t he jet data from HERA [93, |98| and 



Tevatron [9i-[l03. The differences observed between the 
results of A/,(x,Q 2 )/ f t (x,Q 2 ) for KKT12 PDFs to the 
KKT12C PDFs at the input Q\ = 2 GeV 2 and even at 



higher scale Q 2 = 100 GeV 2 , for the u, d and s distribu- 
tions are very small and they are slightly similar at small 
to large x region. 

As Figs.HI^andfTOlshow. the mentioned data sets effect 
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FIG. 6: Our results of KKT12C fit for reduced cross section, & ry Ncfa> Q )? ^ n comparison with HERA combined NC e p reduced 
cross section [IJ] as a function of x for different values of Q 2 bins. The error bars indicate the total experimental uncertainty. 



the gluon and valence quark PDFs at small x, i.e. x < 
10 -3 . The KKT12 gluon and d-valance quark distributions 
both for input and higher scale are dominant at very 
small x, while the KKT12 ratio for u-valance quark PDFs 
are slightly decreases at this region. 



C. Comparison of charm and bottom structure 
functions to the HERA data 

Charm and bottom contributions to the proton struc- 
ture functions, F^x^Q 2 ) and F|(cc, Q 2 ), have been ex- 
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FIG. 7: Our results for reduced cross section, a r ,Nc{x,Q 2 ), in comparison with the recent HI NC ep reduced cross section 
data (6t| as a function of x for different values of Q 2 . 



tracted from Eq. [9] The detailed discussion of corre- 
sponding contributions have been presented in Sec. Ill El 
The NLO (MS) charm contribution F£(x, Q 2 ) to the pro- 
ton structure function in the strict n/ = 3 of the FFNS 
with m c = 1.41 as a function x for different values of Q 2 
and as a function of Q 2 for various x values are shown 
in Figs. [T71 and [TS1 respectively. For comparison we also 
display the charm production data j78l485l | as well. 

Fig. [H shows the NLO (MS) bottom contribution 
F\(x, Q 2 ) to the proton structure function shown as a 
function of x for various Q 2 value in comparison with 



the charm production data [78|, [8Cf . The structure func- 

12, 25, 45 



tion ratio "^,^'^1 as a function of x at Q 2 



GeV 2 in comparison with HI data [82| has been shown 
in Fig. HO] 



D. Comparison to the inclusive jet data 



As indicated in Sec. IIII Al the high-x gluon distribu- 
tion in our global QCD fits is determined by the HERA 
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FIG. 8: The KKT12 parton distributions at input scale Q\ — 
2 GeV 2 as a function of x in the NLO approximation. 




FIG. 9: The KKT12C parton distributions at input scale Q a 
2 GeV 2 as a function of x in the NLO approximation. 



and Tcvatron Run-II inclusive jet data. These data sets 
have the potential to constrain the gluon density in the 
proton when used as inputs to global fits of the pro- 
ton parton distribution functions. In order to illustrate 
the precision of the KKT12C PDFs, we compare the to- 
tal cross sections of some selected jet production pro- 
cesses at the Tevatron Run-II from CDF and D0. Our 
KKT12C NLO QCD predictions for inclusive jet cross sec- 
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FIG. 10: The u-, d-, s-quark and gluon distributions xg at 
Q 2 = 10 GeV 2 as a function of x in the NLO approximation 
in comparison to the results obtained by CT10 0], MSTW08 
1, ABKM10 Ql and GJR08 0J. 
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FIG. 11: The valence quark distributions xu v (x,Q 2 ) and 
xd v (x, Q 2 ) at Q 2 — 100 GeV 2 as a function of x in the NLO 
approximation which have been compared with the results 
obtained by CT10 0], MSTW08 0, ABKM10 [T(J, GJR08 
fl"H and valence analysis of KT08 [26t ] . 
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FIG. 12: Our NLO results for d — u and d/u as a function of x 
in comparison to the results obtained by CT10 @], MSTW08 
@|, ABKM10 [l(J and GJR08 0. The E866 results [HI!!, 
scaled to fixed Q 2 =54 GeV 2 , are shown as the circles with 
statistical and systematic uncertainties. 



tions are compared to the Tevatron Run-II measure- 
ments. Fig. I2T1 sh ows the CDF measured inclusive jet 
cross sections [l 00[ ] . d 2 a / 'dpi^ 1 dy^ et , as a function of p^ 1 
in five different \y-* et \ regions compared to our KKT12C 
NLO pQCD predictions. For presentation, the measure- 
ments in different \y-' et \ regions are scaled by different 
global factors as indicated in parentheses in the figure. 
Overall the KKT12C predictions of NLO QCD gives a good 
description of the data. The D0 Collaboration also re- 
ported their measu rement of the inclusive jet production 
cross section [lOl| . in which a comparison was m ade t o 
our NLO theory ca lculation (with fastNLO code |l04| ) 
based on NL0JET++ [TMEqI] using the KKT12C PDF set. 
The corresponding results of the inclusive jet cross sec- 
tion measurement are displayed in Fig. [22] in six \y^ et \ 
bins as a function of pj? . The comparison shows a good 
agreement between our NLO theoretical predictions and 
the data. 



V. THE LONGITUDINAL STRUCTURE 
FUNCTION F L (x,Q 2 ) 
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FIG. 13: Ratios of Afi(x,Q 2 )/ fi{x,Q 2 ) for KKT11 PDFs 
(fitted to the separate HERA data sets) to the KKT11C PDFs 
(fitted to the HI /ZEUS combined data set) at the input scale, 
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FIG. 14: Ratios of Afjjx, Q 2 )/fi(x, Q 2 ) at the higher scale, 
Q 2 = 100 GeV 2 pi [ml. 



In this section we turn to the perturbative predictions 
for longitudinal structure function Fl(x,Q 2 ). Similarly 
to F^{x, Q 2 ) one can write for Fl(x, Q 2 ) in the common 
MS factorization scheme for n/ = 3 light quark flavour, 

x-^F L (x,Q 2 ) = C L ms® (- qt + ^qt) 

+ 2 - (C LA ® S + C L , g ®g)+ x-'F^ix, Q 2 ) . (15) 

Since the longitudinal structure functions Fl(x,Q 2 ) 
contains rather large heavy flavour contributions in the 
small- x region, a consistent analysis has to be taken into 
account for these effects. The total longitudinal struc- 
ture function can be written as a sum of Fl(x,Q 2 ) = 



F l j* 9ht (x, Q 2 ) + F^ eavy (x, Q 2 ) where heavy quark contri- 
butions sxeF^ avy = Fl(x,Q 2 ) + Fl(x,Q 2 ), and the top 
quark contribution is negligible. 

The extracted parton distributions from present anal- 
ysis have been used to predict Fl(x,Q 2 ). The predic- 
tion for the longitudinal structure functions Fl(x,Q 2 ) is 
mainly determined by the form of the gluon distribution 
extracted from the global QCD analysis. Measuring the 
structure function Fl therefore provides a way of study- 
ing the gluon density and a test of perturbative QCD. 
The study of Fl(x,Q 2 ), led us to extract heavy quarks, 
especially gluon distribution function at low x. 

For completeness, we finally shown our NLO QCD 
standard predictions for Fl(x, Q 2 ) with a rep resentative 
selection of HERA-HI data [gfl M, WH GjJ in Fig. [2J 
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FIG. 15: Ratios of Afi(x, Q 2 )/fi(x, Q 2 ) for KKT12 PDFs to 
the KKT12C PDFs including jet data at the input scale, Q% = 2 
GeV 2 . 
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FIG. 16: The same ratios, A/j(x, Q 2 )/fi(x, Q 2 ), at the 
higher scale, Q 2 = 100 GeV 2 . 



and HERA-ZEUS data [H in Fig. El which we used in 
our QCD analysis. 

The predicted longitudinal structure function 
Fl(x,Q 2 ) as a function of x has been also com- 
pared with the recent HI data [i^] in Fig. [25] Our 
results for F L (x 7 Q 2 ), being gluon dominated in the 
small- x region, are in full agreement with present 
measurements of Fl(x, Q 2 ). 

The longitudinal structure function Fl(x,Q 2 ) as a 
function of Q 2 and for different values of averaged x are 
shown with the most recent HI [69| data in Fig. EH The 
full error bars in this figure include the statistical and 
systematic uncertainties added in quadrature. The inner 
error bars represent statistical error. 
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FIG. 17: The standard NLO (MS) charm quark contribution 
F£(x,Q 2 ) in the strict n f = 3 of the FFNS with m c = 1.41. 
The charm production data are taken from [78l. l80l l8ll. |90| . 
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FIG. 18: The standard NLO (MS) charm contribution 
F£{x,Q 2 ) to the proton structure function shown as a func- 
tion of Q 2 for various x value. For comparison we also display 
the charm production data [78 
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FIG. 19: The standard NLO (MS) bottom contribution 
F%(x,Q 2 ) to the proton structure function shown as a func- 
tion of Q 2 for various x value. For comparison we also display 
the charm production data [78l. l8(i|. 
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FIG. 21: The CDF inclusive jet cross sections measured 
(circles) as a function of Py et for jets with p^f* > 54 GeV/c in 
different \y let \ regions compared to NLO pQCD predictions 
based on the KKT12C PDF (black lines). 
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FIG. 20: The structure function ratio ^i^'^yil as a function 
of x at Q 2 = 12, 25, 45 GeV 2 in comparison with HI data 



FIG. 22: The D0 inclusive jet cross sections in six rapidity 
regions compared to our NLO pQCD predictions based on the 
KKT12C PDF. The data points are multiplied by 40, 20, 8, 4, 
and 2 for the bins \y\ < 0.4, 0.4 < \y\ < 0.8, 0.8 < \y\ < 1.2, 
1.2 < \y\ < 1.6, and 1.6 < \y\ < 2.0, respectively. 
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FIG. 23: Predictions for the longitudinal proton structure 
functio n, F l(x, Q 2 ), at NLO in comparison with HI data [6(il . 
I6ll.l9ll.flil]. 
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FIG. 25: Predictions for the longitudinal proton structure 
function, Fl(x,Q ), at NLO in comparison with recent HI 
data [fiH. 
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FIG. 24: The longitudinal proton structure function, 
Fl(x,Q 2 ), at NLO in comparison with ZEUS data [66| . 
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FIG. 26: The longitudinal structure function Fl{x, Q 2 ) shown 
as a function of Q 2 for different values of averaged x has 
been compared to the most recent HI data for longitudinal 
proton structure function [E^| . The full error bars in the figure 
include the statistical and systematic uncertainties added in 
quadrature. The inner error bars only represent statistical 
error. 
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VII. SUMMARY AND CONCLUSIONS 



The determination of a s is still an open problem in 
PDF fits since it's treatment is closely related to the most 
important processes at hadron colliders such as the LHC. 
This importance is due to the dependance of such pro- 
cesses and their partonic cross section to strong coupling 
constant, a s . We include a s (Qg) as a free parameter in 
our fits, thus the running coupling constant is determined 
in our analysis together with the parton distributions of 
the nucleon; in particular it is closely related to the gluon 
distribution which drives the QCD evolution. Although 
some theoretical groups, like CT10 Q or NNPDF 0,y|, 
fix the a s value close to the updated Particle D ata G roup 
(PDG) average as a s (M|) = 0.1184 ± 0.0007 [ll5| . but 
other theoreticalgroups such as MSTW08 @, ABKM10 
flo| or GJR08 [ll| determine a s as a free parameter in 
the QCD fits. 

The values for the a fl (M|) are obtained in KKT12 fit 
are slightly larger than KKT12C QCD analysis. This is 
due to different data sets which we used in our fits. The 
gluon distributions extracted with and without HERA 
combined data sets play an important role in a s (M|) 
determination. In addition, inclusion of the inclusive jet 
data also effected our a s (M|) values. In our previous 
QCD analysis based on separate and combined data sets, 
which we didn't include the HERA and Tcvatron jet data, 
we obtained a s {M 2 z ) = 0.1185 ± 0.0021 for KKT11 to be 
compared with a s (Ml) = 0.1165 ± 0.0013 for KKT11C at 

the nlo fnlfml . 

Comparison of the present a s (M|) determination at 
NLO with the values obtained by other PDF fitting 
groups has been shown in the Table [TVl The a s (M|) val- 
ues of Ref. [H are slightly larger than the corresponding 
results of both sets of our fits. The KKT12 and KKT12C 
values of a s (M|) are both consistent wit h the world av- 
erage value of a s (M|) = 0.1176 ± 0.0020 [TT^ . 



In the present paper, the most up-to-date data on deep 
inelastic scattering and related processes for the struc- 
ture functions F^x.Q 2 ) with (i = p,d,c,b), Fl{x,Q 2 ), 
xF^{x,Q 2 ) and 'reduced' cross-section data have been 
analyzed in the standard NLO parton model approach of 
the perturbative QCD. We have produced two new PDF 
sets, KKT12 and KKT12C which include up-to-date DIS 
data as well as the HERA and Tevatron Run-II inclusive 
jet cross section data. The KKT12C is obtained using the 
HERA combined data sets while the KKT12 included the 
HERA results as separate data sets. 

Our predictions for all PDFs using QCD fits to the 
data from all available experimental data and the most 
recent HERA combined measurements for KKT12 and 
KKT12C respectively, are in very good agreement with 
the other theoretical models. The strong coupling con- 
stant obtained from our standard NLO analysis in the 
FFN scheme with rif = 3 active light (it, d, s) flavours is 
a s (Mj) = 0.1149 ±0.0019 (KKT12) to be compared with 
a s (M§) = 0.1142 ± 0.0027 (KKT12C). 

A FORTRAN package (grid) containing both our 
sets of standard NLO (MS) parton densities, the light 
F2 9ht (x, Q 2 ) as well as the heavy F^ lea " 1 '(a;, Q 2 ) and lon- 
gitudinal Fl(x,Q 2 ) and xF% structure functions can be 
obtained by electronic mail or can be found di rectl y from 
http : //particles . ipm . ir/links/QCD . htm |l2l| . This 
FORTRAN package covers the following ranges in x and 
Q 2 , 10- 9 < x < 1 and 1 GeV 2 < Q 2 < 10 6 GeV 2 define 
for high Q 2 and/or low x. The QCD analysis based on 
the NNLO approximation is in progress. 



Reference q7(M|) Notes 



KKT12 




0.1149 ± 0.0019 standard approach 


KKT12C 




0.1142 ± 0.0027 standard approach 


MSTW08 


m 


1202 +0 - 0012 
u.izu/ _ 0015 


KT08 


[26] 


0.1149 ± 0.0021 valence analysis 


GJR08 


[39] 


0.1178 ± 0.0021 standard approach 


HI 


[60] 


0.1150 ±0.0017 


CTEQ 


[117] 


0.1170 ±0.0047 


Alekhin 


[US] 


0.1171 ±0.0015 


ZEUS 


[119] 


0.1183 ±0.0028 


BBG 


[120] 


0.1148 0019 va l ence analysis 



TABLE IV: Comparison of the present a s (M|) determination 
at NLO with the values obtained by other PDF fitting groups. 
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